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Аннотация 

Введение. Проектирование систем управления движением мобильных роботов требует построения 
математических моделей. Исследователи неоднократно обращались к этой теме. В частности, опубликованы 
работы, посвященные расчетам мультифизических процессов, моделированию движения различных видов колес 
в тех или иных условиях. Кроме того, рассматривались динамика деформируемых контактирующих тел при 
скольжении, качении и вращении, вопросы автономности и управляемости мобильных роботов. Отметим, что 
динамика и точность позиционирования колесных роботов во многом определяется трением. В литературе не 
представлены исследования динамики робота с дифференциальным приводом с учетом взаимосвязей эффектов 
трения скольжения, верчения и качения на основе теории поликомпонентного трения. Изыскания в этой области 
способны выявить новые динамические эффекты. Основываясь на полученных таким образом данных, можно 
заняться повышением точности позиционирования при построении математических моделей управления. Цель 
представленной работы — исследовать движение автоматического устройства с дифференциальным приводом с 
учетом трех моделей контакта: неголономного, кулоновского трения и поликомпонентного. 

Материалы и методы. В качестве базовой приняли схему двухколесного робота с дифференциальным 
приводом и безотрывным движением по опорной поверхности. Перемещение устройства обеспечивалось 
программным управлением. Динамику описали в виде уравнений Аппеля. Для вычислений использовали 
математические модели, по-разному учитывающие трение. Координация действий механизма формировалась на 
динамическом уровне. Управляющие воздействия — моменты двигателей колес. При визуализации исследуемых 
моделей задействовали встроенные численные методы системы \оЙат Мафетайса («Вольфрам математика») 
с минимальной точностью 10-5. 

Результаты исследования. При построении математической модели определили равенства для угловых 
скоростей колес. Учли наличие площадки контакта и вывели уравнения динамики робота с дифференциальным 
приводом. Элементами системы были проекции силы и момент, показатели верчения платформы, массы, угловые 
ускорения и инерция колес. Показали, как формируются управляющие воздействия в рамках неголономной 
механики. Описали модель двигателей, которые создают момент управления на ведущих колесах. Вывели 
решение как взаимосвязь индуктивности проводников электродвигателей и работы электропитания. Детально 
рассмотрели три модели, описывающие динамику робота с дифференциальным приводом. Первая — 
неголономная. Вторая и третья включали систему уравнений динамики робота с дифференциальным приводом 
для общего случая, при наличии площадки контакта. При этом во второй модели игнорируется время 
переключений в двигателе и задействуется трение Кулона. В третьей ввели параметр, определяющий скорость 
переходных процессов в двигателе, и задействовали разложения Паде. Это модель с поликомпонентным трением. 
Итоги расчетов показали в виде графиков. На них исследуемые модели визуализировали в виде кривых разного 
цвета. Сопоставление графиков показало, в каких случаях после завершения переходных процессов управление 
обеспечит требуемую точность. 

Это модели 1 и 2. В модели 3 программное управление генерирует ошибку в угловой скорости вращения 
платформы. Данную ошибку нельзя спрогнозировать в рамках 1-Й и 2-й модели. Во всех рассмотренных системах 
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скорость скольжения колес в поперечном направлении падает до нуля. Получено и обосновано условие 
безотрывного движения опорного колеса. 

Обсуждение и заключение. Программное управление допустимо в моделях, не учитывающих трение колес при 
одновременном скольжении, верчении и качении (общий случай пространственного движения). Однако важно 
учитывать взаимосвязь данных процессов и поликомпонентного трения. Это необходимо для более точного 
выполнения роботом программных движений. Установлено, что программное управление в модели, 
учитывающей трение верчения и качения, приводит к отклонениям от программных значений угловой скорости 
платформы. Полученные результаты можно использовать при построении системы управления 
с прогнозирующими моделями. 


Ключевые слова: программное управление роботом, управление в неголономной механике, модель трения 
Кулона, поликомпонентное трение, ошибка угловой скорости вращения 
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Введение. В современном мире растет востребованность автоматических мобильных колесных аппаратов, 
активно исследуются возможности координации их работы. Значительный научный и прикладной интерес имеют 
вопросы автономности и управляемости мобильных роботов. В частности, рассматриваются подходы к расчетам 
мультифизических процессов, моделированию движения конкретных видов колес в разных условиях. Некоторые 
изыскания фокусируются на динамике деформируемых контактирующих тел при скольжении, качении и вращении. 
Результаты этих исследований находят применение в робототехнике [1-3]. Управление мобильными колесными 
аппаратами рассматривается на примере электрических самокатов «Сигвей» (Зез\мау) [4]. 

Мобильные механизмы с простыми и эффективными системами управления движением колес нужны 
в первую очередь для работы в узких пространствах [5], для проведения исследований в сложных 
промышленных условиях [6]. 

В [711] показаны возможности точного моделирования динамики колесных аппаратов с учетом скольжения, 
верчения и качения колес. 

При этом в литературе не описаны явления, возникающие при комбинации эффектов трения. Цель 
представленной работы — исследовать влияние поликомпонентного трения на динамику робота 
с дифференциальным приводом. Задачи: построить математическую модель движения робота 


с поликомпонентным трением и смоделировать программные движения. 
Материалы и методы. Рассматривается схема двухколесного робота с дифференциальным 
приводом (рис. 1). Предполагается безотрывное движение колес радиуса К по опорной поверхности. 
Хх] 
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Рис. 1. Кинематическая схема робота: а — платформа; 6 — одно из ведущих колес 


Салимов М.С. и др. Программное управление движением робота с дифференциальным приводом при разных моделях трения 


Центр тяжести платформы С смещен относительно геометрического центра платформы А по оси Х на 
расстояние Й, ‚ по оси У — на расстояние й,. Точки С; и С, — центры масс колес (индекс 1 соответствует 


левому колесу, индекс 2 — правому). 

Вводятся следующие обозначения: т —щ масса платформы робота; с, — момент инерции платформы 
относительно центра масс С; Ох, у — неподвижная система координат; х,, у, определяют положение точки 
Ав О,х,у,; у — угол курса платформы; ф, — углы поворота колес (индекс 1 соответствует левому колесу, 
индекс 2 — правому); АС =АС, =Ё — расстояние между колесами; М,,М, — моменты двигателей, 
приложенные к колесам. 

Предполагается, что на опорное колесо 4 (рис. 1) действует только нормальная реакция опоры. Во время 
движения опорное колесо располагается на оси Х на расстоянии А, от центра платформы (А, < Г). 

Результаты исследования 

Математическая модель динамики робота. При описании движения робота с дифференциальным приводом 


в рамках неголономной механики применяются выражения, которые учитывают ограничения, связанные 
с отсутствием скольжения колес по отношению к поверхности: 


0х, =00х, =0,ох,, =0, (1) 


Ок 


1х 
где Эк Ок,х> Оку ОК, —Щ_ проекции на оси Хи? скоростеи точек К, иК, › являющихся центрами площадок 
контакта колес с опорной поверхностью. 

Отметим, ЧТО ДЛЯ неголономной модели предполагается точечный контакт колес с поверхностью, 


т.е. площадка контакта каждого колеса вырождается в единственную точку контакта. 
Рассмотрим выражения: 


ок. =0-©ОЁ- Ав, 
ок =о+0- К®,, (2) 
Ок б Оку =": 


С их помощью можно получить равенства для угловых скоростей колес: 
о- ОР 5+ ор 
в, = ‚в, = : (3) 
К К 
В общем случае, когда предполагается наличие площадки контакта, получим уравнения динамики робота с 
дифференциальным приводом: 


т, [6-0-9(и +1,0) |+2т, (5-Ои) = А, соза, + А, соза, — Е эт а, — Е, та, 
т, [#+л.0+9(5-л,9) | +2т, (й+05) = Е, эта, + Е, та, — Е. 0$ а — Р,.с08 а,, 


(с. +2. +т, (+ )+2т,Р )+т, [(и+95)л, —й,(6-©и)] = 


4 
Мс +М.. +(5, со5а, — Е зш а, —Ё, с0за +Ё 5 о, ) Е, _. 
Ль®, =М,-М„ -(Е, соза, - Е, зто, К, 
Люб, =М,-М„, (=, сова, -Ё\, зто, ) К. 
Здесь К,Ё, — проекции силы трения на направления скоростей скольжения центров пятен контакта колес; 


Е, 5, —_ Проекции силы трения на направления, перпендикулярные скоростям центров пятен контакта колес; 


М.,М.,— слагаемые момента трения при верчении для двух колес; @ — угловое ускорение верчения 
платформы; &,,6, — угловые ускорения колес, проецируемые на их ось вращения; т, — масса каждого колеса; 
Ук, и Лк. — моменты инерции колеса относительно соответствующих осей. 


Для формирования управляющих воздействий будем использовать модель динамики робота, описанного в 


рамках неголономной механики [12]: 
ь 2Л М+М 
т [5-9 ол,о |+2т =, 
| | К К 
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т, Г(+^.5 +9(5-л,9))л, —^, (-л,0-0^,0) | у 


5) 
р. (М,-М,)г ( 
2 2 2 — 2 1 
(ле. +2. +т (в +®)+2т,Р) +25 кто —— 
Ее также можно представить в виде: 
м, +М м,-м,)Е 
Е = 1 те 2 „Мо — ( 2 ы , (6) 


где Ё, и М. — левые части уравнения (5). 
Модель двигателей на ведущих колесах. Опишем модель двигателей, способных генерировать момент 
управления на ведущих колесах М,,М,. 


Уравнения движения колесных двигателей постоянного тока можно выразить формулами [13]: 


4(1.) 


в в. @) 

М =сж, (^ =) 
Здесь Г, —_ индуктивность проводников электродвигателей; в — сопротивление электрическому току в цепи 
двигателя; и, —_ электропитание для двигателя к 1, —щ ток, проходящий через цепь якоря двигателя; СС —= 


конструкционные постоянные двигателей. 
Предполагается, что двигатели ведущих колес имеют одинаковую характеристику. Из (7) видно, ЧТО 


ре. (8) 


Учитывая (8), в первом уравнении из системы (7) получим дифференциальное уравнение относительно 
момента М, : 
Е... В 
е ыы 
оп о. —с0,. (9) 
2 


ё 


К 
Уравнение (9) разделим на множитель — и получим: 


с2 
Р.С. с сс 
- -М,+М, =-—® 0, - 2 ,. (10) 
[6 К К 
Введем обозначения: 
Ес сс [6 
ПЕК =, К, =. (11) 
с, К у У: 
В таком случае управляющие моменты можно описать: 

ТМ, +М, =К.Ц, -К\о,, (12) 
где К,К, — константы; 1, — параметр, определяющий скорость переходных процессов в электрической части 
двигателей. 

При игнорировании времени переключений в электромеханической системе двигателя: 


М, =К,И, -Ко,. (13) 


Авторы [14] получили похожие формулы. 

Исследуемые модели. Рассмотрим модели, описывающие динамику робота с дифференциальным приводом. 

Модель 1 описывается уравнениями (6) и (12), когда длительность переходных процессов полагается 
исчезающе малой 1 =0: 


: 2/7, Ки, -К КИ, -К 
т [6-л,о-Ол,О ]+2т,5+ Ку о = КОД 10, + 50. > 
1 р; р: к В? В 


2 


т, [(+^.0 +9(5-л,о))л,—й, (5-1,0 -9^,9). +2. 0+ 


14 
+0 (о. +2.Л. +т, (2 +й, ) + 2т, 1?) = (к, ВВ р : _ 
К 
Здесь ®,,®, вычисляются по формулам (3). 
Для описания модели 2 применяются уравнения (4), (13), а также модель трения Кулона: 
Е =-ЛМ,Е, =0,Мс =0. (15) 
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Модель 3 представлена уравнениями (4), (12). Кроме того, учитывается поликомпонентное трение, 
представленное в [12] с помощью разложений Паде: 


21 -К. 1, со [ ма с) 
оо, ЕЛЕЫ - со-[ “) 


Е =-ЛМ№, ° = › (16) 
21 -К. 1, сиз[ А-а) 
А Аа 
0 
Кпо 
Е, =- ЛМ, з (17) 
21 —К, 1, со п-ов 
п 
КГ, т. сон с) ок, Нок, +] 
п 
ешь [ А-а) 
' 27 
2 2, со-[ с) ох, +0 
рт Га 
Ме =-2 1, . (18) 
в - 
А, 08| ——а | |2. +9х, +] 
Ь 
М, = ль". (19) 
| || 


Здесь [,1[,,/, — постоянные величины, которые находятся по формулам в [12]; и — угол, характеризующий 
направление линейных скоростей скольжения относительно точек контакта колес КК; К — 


7 


коэффициент, знак которого зависит от направления качения. 
Для определения программного управления следует указать законы движения робота 5“ (1),О“ (1). По (5) 


можно выразить величины Ё, = у (0, Е = т (г), Мо=М ь (г) в зависимости от времени. Используем 
известные значения Р., М. для вычисления управляющих напряжений (/ и (0, , подаваемых на приводы колес: 


В М К М. 
Ч, ре Пка ко Це = пк -ь Ко ь (20) 
- К (2 7Й К, \ 2 р 


м м 
Движение системы с заданными начальными условиями 0(0)=0,5—,и(0)=0,2—,©(0) = 0,2 
с с с 


смоделировали в математическом пакете «Вольфрам математика» (\МоШат Мафетайса). 


т 1 
Рассмотрим 5“(0=0,5 —, 0“ (1) =0,3—. Для интеграции системы автоматически выбирали шаг 
[6 [6 


интегрирования и задействовали встроенные численные методы У\УоШат Мафетайса с точностью результатов 
не менее 105. 


У(В, м/с 


Рис. 2. Зависимость продольной скорости у. Зеленая кривая — неголономная модель 1, 
красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением. 
Синяя и красная кривые совпадают 
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и(1), м/с 
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Рис. 3. Зависимость поперечной скорости и. Зеленая кривая — неголономная модель 1, 
красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением. 
Зеленая, синяя и красная кривые совпадают 


©(0, 1 
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Рис. 4. Зависимость угловой скорости робота. Зеленая кривая — неголономная модель 1, 


красная кривая — модель 2с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением 


Из рис. 2 и 4 видно, что управление гарантирует точную реализацию программных законов движения после 
завершения переходных процессов в моделях | (неголономная) и2 (с трением Кулона). В модели 3 
(поликомпонентное трение скольжения, качения и верчения) после завершения переходных процессов из-за 
программного управления возникает постоянная ошибка в угловой скорости вращения платформы. Во всех 
рассмотренных моделях скорость скольжения колес в поперечном направлении быстро снижается до нуля. 
Случаи совпадения зависимостей, полученных для разных моделей, показаны на рис. 2 и3. 

Для проверки выполнения условий безотрывного движения опорного колеса запишем теорему изменения 
кинетического момента системы в проекциях на ось 7: 


Лк, (в, +6.) =М,+М,-М, -М, -М, В, +т, эй, -(Е 0$ 0, — 2 зто +, соза, -Ё, уто, ) К. 


ьЬ 
Согласно (4): 
М+М, =-М, А, +т,8й,. 
Из последнего уравнения получаем условие безотрывного движения опорного колеса: 
тай, —-(М,-М,) > 0. 
С учетом (5): 


а. , 
т)0+ Вт [и,Ф+ 1,0] >0. 
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Ускорения робота не превышали 1 м/с?, поэтому данное неравенство выполняется в каждый момент времени. 

Обсуждение и заключение. Итоги моделирования позволили сделать ряд выводов. Программное управление, 
основанное на неголономной модели 1, приемлемо для роботов, конструкция которых не учитывает контактное 
трение колес, возникающее при комбинации скольжения, верчения и качения. Однако при формировании 
управляющих воздействий для достижения более точного следования программным законам движения 
необходимо учитывать модель поликомпонентного трения. Полученное решение можно использовать при 
построении системы управления с прогнозирующей моделью (то4е! рге@сйуе сопёго|. Такая разработка будет 
предметом дальнейшего исследования. 
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